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fvq Abstract. In this work we study a family of bosonic lattice models that combine an 

on-site repulsion term with a nearest-neighbor pairing term, ^YliU i) ^I^I + H.c. Like the 
(-H original Bose-Hubbard model, the nearest-neighbor term is responsible for the mobility 

^ I of bosons and it competes with the local interaction, inducing two-mode squeezing. 

However, unlike a trivial hopping, the counter-rotating terms form pairing cannot be 
studied with a simple mean-field theory and does not present a quantum phase transition 
r^ in phase space. Instead, we show that there is a cross-over from a pure insulator to 

O^ long-range correlations that start up as soon as the two-mode squeezing is switched on. 

We also show how this model can be naturally implemented using coupled microwave 
resonators and superconducting qubits. 
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1. Introduction 

The physics of lattice bosons is an interesting topic that has experienced several revivals 
associated to the proposal and realization of new experimental setups. These include, for 
instance, the study of phase transitions in Josephson-junction arrays [H |2], the ground- 
breaking studies of quantum phase transitions with bosonic atoms in optical lattices [3] 
and, more recently, the proposals for polariton physics in coupled cavities and photonic 
systems [H El E] . The previous represents a curious round trip where a very interesting 
physical model is introduced in the context of Condensed Matter Physics, we learn how 
to tame it and control it with the tools of Quantum Simulation and the final and most 
promising platform for studying its physics seems to be, once more, a superconducting 
setup [HE]. 

The use of superconducting quantum circuit setups to implement bosonic lattice 
models has various advantages over atomic implementations, but the most important one 
is the access to a wider variety of interactions. Using photons as particles we can easily 
engineer nonlinearities [H El E] , but we also have access to controlled dissipation [9l [10] , 
easily customizable geometries [S] and effective dispersion relations [H] and, most 
important, the possibility of increasing the coupling to a point in which counterrotating 
interactions become relevant [121 [13] • 

In this work we explore the possibility of engineering exotic interactions in circuit- 
QED, studying models with photons that do not preserve the number of particles. Our 
study is centered around a model of coupled cavities with pairing and on-site interaction 
{h = l) 

if = ^ ^ (e''^'^a\a] + H.c.) + f^ ^^ a\a\aiai + ^ Uia\ai. (1) 

In contrast with the original Bose-Hubbard model 

H = -^-Y^ (e'^'^alaj + H.c.) +UJ2 44<^^ciu (2) 

{hi) « 

our problem facilitates the mobility of bosons through a pairing term, a\ap that does not 
preserve the number of particles. We will show that the competition between pairing, g, 
and on-site repulsion, [/, does not have a phase transition associated to it. Instead, for any 
small value of g we find a cross-over mechanism that establishes long-range entanglement 
and squeezing in the lattice, through a process that cannot be described by a trivial 
mean-field theory. 

Our work is structured in three sections. The first two sections study this model in 
two different regimes. In Sect. [2] we adopt the regime f/ = of non- interacting bosons and 
derive the squeezing and entanglement properties of the resulting array of linear cavities. 
We will show that the ground state in momentum space is a collection of paired photons 
with different momenta. The paired momenta can be controlled using the phase 0jj in 
Eq. ([T]) and the whole system approaches criticality when g — to. The competition between 
this long-range entanglement and the local interaction is the topic of Sect. [3j There we 
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will add a nonzero value of U and study the problem using the methods of infinite Matrix 
Product States or iTEBD [13]. We will show that now U prevents a breakdown of the 
model for large g and that it also suppresses entanglement. However, unlike in the Bose- 
Hubbard model, the transition from the insulator regime to the squeezed multimode 
system does not happen through a phase transition but through a cross-over. In Sect. |4] 
we propose a quantum simulation of Eq. ([I| using an array of superconducting microwave 
resonators which are coupled through periodically driven SQUIDs [15J. We will consider 
important issues, such as state preparation and tomography of the correlations that were 
estimated in previous sections. Finally, in Sect. [5] we summarise our results and discuss 
possible extensions of this work. 

2. Linear Cavities 

In this section we will study a uniform model of coupled cavities {ui = u) without on-site 
interactions (f/ = 0) 

Just as in the Bose-Hubbard model, our aim is to find the state of the bosons in the limit 
in which they are free particles. Since this is a quadratic problem in Fock operators, we 
will use the formalism of Gaussian states, analyzing the entanglement that arises in the 
system which, as we will see, is not condensed but forms a multimode squeezed state. 



In Section 2.1 we will first Fourier transform Hn in order to write it in momentum 



space, where we will see that it can be expressed as the direct sum of two-mode 



Hamiltonian. Then, in Section 2.2 we will be able to diagonalize it and study its stability 



conditions. In Section 2.3 we will analyze the effect of adding a phase to each site 



that grows linearly with the lattice. In Section 2.4 we will study entanglement between 



photonic modes, i.e. between modes in momentum space. 

2.1. The reciprocal lattice 

For simplicity we will assume that the model inEq. (pi) describes a one-dimensional 
lattice with periodic boundary conditions, and assume (pij = 0. Due to the translational 
invariance of the problem, it makes sense to rewrite the same model in Fourier space. We 
perform the transformation 

«^ = ^E^"""V- (4) 

^ ^^ k£i 

Here the lattice sites run over the indices i = 0, . . . ,N — 1, and qk = 2TTk/N G 
(— vr, vr] is the quasimomentum and we choose the quasimomentum index to run over 
/ = {- [f - 1] , . . . , [f ] } if AT is even or over / = {- [^] ,..., [^] } if it is odd, 
forming a symmetric Brillouin zone. This unitary transformation transforms Hq to a sum 
of Hamiltonians, Hq = ^^ Hq with 

H^o'^^ = uj{ng + n^q) + g cos(g) {bgb^q + H.c.) . (5) 
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Note how, except for Qk = and q^ = ii (if A^ is even), H^ represents the usual two- 
mode squeezing interaction between bosonic modes. In the following we will analyze the 
entanglement that results from this coupling, neglecting the point-like singularities at 
Qk = 0, TT, which do not add much to the physics. 

2.2. Two-mode squeezing 

We diagonalize the momentum Hamiltonian dsl) using the squeezing transformation 



I {c,b,b., - e.&;&i,) 



e, = r,e'^' (6) 



Sq{Q = exp 

obtaining an operator 

^i^) = Sl{QHi;^^S,{Q = e, {b% + bl^b^,) +e,-u. (7) 

The dispersion relation of this problem is 

Eg = a/w2 - g2cos{qy. (8) 

In order to warrant the stability of the setup we must constrain the photon-photon 
interaction to be \g\ < u, because otherwise there would be no vacuum for this problem. 
Doing so, the squeezing parameter is defined as the oscillating function 

tanhr,= '^^"^^^^' (9) 

CO 

which is maximum around g = and ivr. Since we have chosen 0jj = 0, the parameter 
^q is real, and the phase (pg just absorbs its sign. 

2.3. Momentum boosts 

In the full model we considered the possibility of having a site-to-site phase difference, 
(pij = (pi- If this phase is constant, it may be eliminated through a global phase 
transformation of the bosonic operators. Instead we will consider it to grow linearly 
through the lattice. Since we wish to preserve translational invariance and periodic 
boundary conditions, we find that the dependence of the phase must be 

0, = 7r^(2z + l), rf,zG{0,...,iV-l}. (10) 

We can engineer a unitary gauge transformation that absorbs this phase 

a, = e^-'^^/^5,. (11) 

The original Hamiltonian Hq is expressed now in terms of {Sj, a-}, but with (pi = 0. The 
difference is that the quasimomentum operators for this transformed Hamiltonian 

correspond to the quasimomenta q^ — nd/N from the original model. In other words, the 
phase (pi = 2iTd/N corresponds to a momentum boost of amount —ird/N. Note that to 
achieve a similar effect in the Bose-Hubbard model we needed a Peierls transformation 
with uniform phase. 
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Figure 1. (a) Logarithmic negagtivity Ej^ as a function oi g/u! and q for 0^ = 0. Every 
quasi momenta q^ is entangled with q^ = —q-. The logarithmic negativity diverges at 



(b) Similar plot for 0, = 7r^(2i + 1) (see Eq. lOl, with 



q = and q — n ior g -^ 

d = 0.2N. Every quasi momenta q- is entangled with g_)_ = — (/_ — 2nd/N, and the whole 

Brillouin zone is boosted by —nd/N. 



2.4- Entangled Photon Pairs 

Since the Hamiltonian splits into a sum of commuting terms, each of them acting on a 
separate pair of momenta, we may regard the total density matrix as a tensor product 
P = ®Pq of Gaussian states, pg, for each of the momenta. Each of these density matrices 
is best characterized by the covariance matrix (CM). 

The CM is written in terms of first and second order statistical moments using 
the field quadrature operators Xq = {bq + 6j)/-\/2 and Pq = i(6j — bq)/\/2. Defining 
Qig) _ ^rf.^^ p^^ 2;_q, P-q)^, the CM is the 4x4 matrix that results from the expectation 
values Cnm = \{QriQm + QmQn) — (Qn) (Qm) ■ In Complete agreement with any quantum 
characterization of two-mode Gaussian states [T6], we find a structure of 2 x 2 blocks 




C 



where a = /3 = | cosh r^ ■ I2 and 7 
correlations between modes. 



(13) 



|sinhrq. Out of these terms, only 7 encodes 



According to the Peres-Horodecki-Simon criterion [T71 UHl [19] the positivity of the 
partial transpose C^^ = A'^CA, A = diag (1, 1, 1, — 1) represents a sufficient condition 
for separability of the two-mode state. The entanglement monotone associated to this 
criterion is the logarithmic negativity 



(14) 



I 4 

En = --^log2[min{l,2|z/i|}] 
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where z/, G {— z^±, +i^±} are the four symplectic eigenvalues of the covariance matrix 



,'Ai±v/A?-4A2 1 litanhr, 
'^ = \J 2 = 2V lTtanhr/ ^^'^ 

where we substituted the invariants Ai = det a + det /3 — 2 det 7 and A2 = det C for their 
value in this particular problem. 

The immediate consequence is that the negativity is nonzero for any amount of two- 
mode squeezing 

P l, A + l^cos(g)/w|\ 

EN,q = - log2 . ^^^ . (16) 

2 \l-\gcos{q)/uj\J 

In other words, as it was intuitively expected the state is not separable and the 

entanglement grows with increasing coupling strength, g. 

In Fig. [T]we show a typical distribution of the entanglement on the first Brillouin zone, 

both without and with the phase 0j. In both cases we have pair entanglement between 

different momenta, which is maximal around to points. The maxima are displaced with 

the help of the momentum boost (pi. Using 0j we can switch from entangling photons 

that travel in opposite directions but carry the same energy, to entangling co-propagating 

photons that differ in frequency, as is the case of 0j 7^ in Fig. [1]d. Note also the similarity 

between the shape of the entanglement and that of the squeezing (c. f. Eq. ([9])). 

3. Pair-Hopping Tunneling 

So far we have seen that in the linear case, the counterrotating terms create entanglement 
and long range correlations. We want to analyze how these correlations compete with an 
on-site interaction such as the U term from our full model in Eq. ([I]). 

We already know from the physics of the Bose-Hubbard model (cf. Eq.(|2])) a similar 
competition |2T]: a kinetic term which preserves the number of particles fights for their 
delocalization against an on-site repulsion that hinders mobility. The result is a quantum 
phase transition between a superfluid regime, for t ^ U, and an insulator regime, U ^ t, 
which has been repeatedly observed using bosonic atoms in optical lattices [3] . 

In the model from Eq. ([I]) the kinetic term is implemented with a pairing term that, 
by allowing the creation and destruction of pairs of particles in nearby sites, it also allows 
the mobility of particles. This is obvious from the previous results, the energy bands and 
the delocalized wavefunctions that result. It is obvious that such a term would compete 
with the on-site repulsion U, but it is not clear on which terms this competition happens. 

A trivial inspection of Eq. (fTl) reveals that, should we try a similar mean-field 
treatment as for the Bose-Hubbard, that is replacing a with a nonzero expectation value, 
(a), we would obtain the same phase diagram as for the original model. More precisely, 
assuming a product state ansatz \^p) = ®^o^ \^)^ ^^^ energy per site reads 

^(^).=o = f [(a)^ + H.c.]+f/|(a)r, (17) 

which is minimized for a nonzero expectation value of a. Yet we know from the linear 
model that (a) = always because of the creation and destruction of pairs of particles. 
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(a) , . (b) (c) 




Figure 2. Phase diagram of the free energy H — /iN, for the model in Eq. (fTl). We 
plot the (a) number of photons, (b) fluctuations of the photons and (c) the nearest- 
neighbor pairing term, vs. the chemical potential fi and the photon-photon pairing g. 
The usual boson delocalization {a\aj) is zero over the whole diagram. An occupation 
number cut-off n„i,ax = 10 for the states {ik\ (cf. Eq. ([l8|) and x = 20 has been choosen. 



This opens the door to a different type of phase transition, but it also points out the need 
of a more rigorous and more sophisticated treatment of the problem. 

In the following pages we will show that indeed there is no such phase transition, but 
instead we find a continuous growth of entanglement along the lattice as the pairing g 
increases and U decreases. This is proven using two different tools: the iTEBD method 
developed in Ref. [22] and further explained in Ref. [11], and perturbation theory in the 
weak coupling regime, g <^U . 



3.1. Bosonic phase diagram 

The iTEBD method [22j poses a translationally invariant tensor-product state ansatz for 
a one-dimensional quantum system. Such an ansatz may be qualitatively written as the 
contraction of four tensors, which are repeated with period two along the lattice, 

|^) = ^...r^'^«r^'^T^'^^..|...zo,^i,^2,...). (18) 

Here, the [ik] represent the physical states of the fc-th lattice sites, which in our case is the 
occupation number of the cavities, cut-off to a reasonable number: ik = 0,1,2... Umax- 



The tensors F have dimension x^ 



Ur, 



+ 1), where x is another cut-off, this time for the 



amount of entanglement that the ansatz may host. 

Taking into account the obvious limitations of the ansatz, it is nevertheless possible to 
get qualitative and even quantitatively accurate answers to many of the physical properties 
of the ground state. We have applied the ansatz to the study of the free energy function 

F = H^=o - /iiV, (19) 
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Figure 3. Perturbation theory results (dashed) and iTEBD expectation values (solid) 
for the full model in Eq. (fTl), close to the insulating phase {g/U <C 1). We plot the (a) 
number of photons, (b) fluctuations of the photons and (c) the nearest-neighbor pairing 
term, vs. the photon-photon pairing g. The cut-off n„iax a-nd x (cf- Eq. (18l) are the 
same as in Fig. [2] 



where /x is a chemical potential that favours the injection of photons in the model ([I| and 
N is the total number of photons. The choice of the free energy instead of the Hamiltonian 
is based on the need of a comparison with the Bose-Hubbard model. 

Minimising the energy (F) over the set of iTEBD states with fixed bond dimension, 
X, we have produced the pictures shown in Fig. [2| where we plot the density, the number 
fluctuations (An)^ = ((ajaj)^) — {alai)'^ and the pairing between nearest neighbors, all 
for the ground state. In these pictures we appreciate that the lobes of the Bose-Hubbard 
model are present in the growth of the number fluctuations [cf. Figs. ^]. However, the 
system no longer exhibits plateaus of the density, and there is no sharp phase transition 
between the insulator and the squeezed or entangled regime. The order parameter of a 
typical superfluid, which was the largest eigenvalue of the single-particle density matrix, is 
now zero because {alaj) = everywhere. Instead, the new order parameter is associated 
to the squeezing, which we identify with a very simple observable, the pairing, {aloj). 
This correlator increases smoothly with the interaction g/U, as shown in Fig. ^. 

The continuity of the transition can be further confirmed by studying weak 
perturbations around an insulator state. We start with a bare sate that contains exactly n 
photons per cavity |u°) = (3)j=o \^)iy and apply perturbation theory up to second order, 
obtaining (cf. Sec. Appendix A[) 



(alai) 



n + ^ (/„%! - fn) + O {g') 



n' + ^- [i2n + !)/„%, - (2n - l)/^^] + O {g') 



\ ili U^i 



,9 



nfn - {n + 1)/, 



n+lj 



0{9') 



{a -a. 



i+l 



) =0 + O{g') 



(20) 

(21) 
(22) 
(23) 



Bose-Hubbard models with photon pairing in circuit-QED 
(a) , , (b) 



0.1 0.2 0.3 0.4 



An 

0.1 0.2 0.3 0.4 0.5 



(c) 

0.05 0.1 0.15 0.2 




Figure 4. Phase diagram of the driven systeni H + ^^ r2(aj + a\), for the model in 



Eq. (24). We plot the (a) number of photons, (b) fluctuations of the photons and (c) 
the nearest-neighbor pairing term, vs. the chemical potential /i and the photon-photon 
pairing g, taking as units uji = U = 1. The usual boson delocalization {a\aj) is zero over 
the whole diagram. 



expressed in terms of the function fn = n [2f/(n — 1) — /i]~^. 

Using these expressions we can approximate all observables, including the density, 
the number fluctuation, etc., and compare with the exact results. This is done in Fig. [3} 
where we find an excellent agreement all the way up to (7 = O.lf/, for a variety of filling 
factors. It is important to remark that perturbation theory confirms the lack of superfiuid 
order, accompanied by the fast growth of the nearest-neighbor squeezing (the quantity 
that grows fastest of all for small g). 

3.2. Photonic phase diagram 

The previous approach relied on a Condensed Matter Physics phase diagram that 
combines the interaction with the chemical potential. This approach is needed when 
studying real particles because the number of atoms and electrons in low-energy 
experiments is a superselection rule. However, for the setups that we have in mind, which 
implement bosons using photonic or polaritonic excitations, there is no such conservation 
rule. Instead, the equilibrium state of the system is usually a vacuum, squeezed or not. 

It is for this reason that, in order to recreate an actual experiment, we have to 
simulate a different type of problem: one with a mechanism that actually injects particles 
in the experiment. For this purpose we implement a very simple term 

i^drivcn = H + Y, ^(«* + «1), (24) 



that consists on an external displacement of the cavity. This displacement, of strength 
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Figure 5. iTEBD expectation values for tlie number of photons per lattice site. We 
plot the number of photons for (a) the driven model in Eq. p4| ), vs. the photon-photon 
pairing g and the displacement il, (b) the full model in Eq. (nj), vs. the photon-photon 



pairing g and the chemical potential /i. The cut-off Umax and x (cf- Eq. (18l) are the 
same as in Fig. [2] 



Vt/ ^/uj in the non-interacting case, is implemented through an external source, which in 
our case can be a time-independent flux source (i.e. an inductor). 

An unfortunate consequence of this formulation is that we lose the lobe structure. 
For fixed Wj = 1, f/ = 1, Fig. |4] shows the same observables as before. We appreciate that 
the only insulating region is now around the vacuum, aX g = VL = Q. Pairing continues to 
be a good order parameter, but it has to be properly defined, eliminating the equilibrium 
values due to the injection, (ajaj^_^) — {a\){a\^]). If we do so, then the pairing is relatively 
insensitive to Vt and the average number of photons, and grows rapidly with the coupling 
strength g. 



4. Quantum Simulation 

The model that we have studied so far ([I| is not likely to appear in nature: the 
counterrotating terms normally appear together with rotating terms, in a way that these 
ones are dominating the dynamics. However, the fact that certain interactions are not 
commonly available does not mean that we cannot study them in a physical setup. Instead, 
we can find a highly tunable quantum mechanical system and implement such interactions 
by controlling its dynamics, in what is known as a "quantum simulation" [See Ref. [23] 
and accompanying articles]. In essence, this is the spirit behind all recent proposals about 
coupled cavities and polariton-type physics [H El E], and it will be the way we suggest to 
implement those models. 

For the simulation of our coupled- cavity model we suggest using a one-dimensional 
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(a) 



(b) 




Figure 6. (a) One dimensional array of resonators coupled by means of a 
superconducting ring coupler, (b) One dimensional array of circular resonators coupled 
by dc-SQUIDs. Both lattices are folded in a ring in order to keep periodic boundary 
conditions. Inductive coupling by Josephson Junctions and intersected loops allow to 
engineer bipartite interaction between nearest neighbors that is governed by counter- 
rotating terms. 



setup of coupled microwave resonators. As shown in Fig. [6} such setup consists basically 
on a number of superconducting segments (coplanar waveguides or striplines) joined by 
a nonlinear element called SQUID. Each piece of cable supports a number of quantized 
standing waves, of which we will only focus on the fundamental mode, ujia\ai, with a 
frequency that will change from cavity to cavity. As shown elsewhere |15j, by joining the 
cables through SQUIDs, we actually implement a tunable coupling element 

H = ^uJiajai + ^g[^i{t)]{al + ai){al^-^ + a^+i), (25) 

i i 

which can be controlled using the magnetic flux $ that runs through the Josephson 
junction loop. 

To activate the counter-rotating (rotating) terms in the previous equation, we only 
need to drive the flux on each cavity periodically [T3] 

^i{t) c:iai + hi cos{uit + ipi) (26) 

with arbitrary a^, hi and phase. If the magnetic flux $ is driven with z/j being on resonance 
with the counter-rotating terms, the propagation of the rotating terms may be averaged 
out and vice versa. For this we only need to fabricate cavities with alternating frequencies 

cj. = cjo + (-l)'(5a;/2, (27) 

and choose a driving vi = 2ujq (alternatively chose i^i = 6u){—iy in order to activate 
rotating terms). 

Note also that in addition to the resonance, we can also control the phase. The 
time-origin of the driving propagates through a unitary transformation to the effective 
model that results, creating the full pairing term from Eq. ([I]). 

The last two ingredients in our system are the interaction and the possible external 
driving. The last one is easy to achieve, because the inductive or capacitive coupling 
between a cavity and a nearby cable will introduce the possibility of driving photons in 
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and out of the system. The nonhnearity is a bit trickier, as it involves a nontrivial photon- 
photon interaction inside the cavity. We expect that such terms could be ported from 
existing proposals based on nonlinear resonators [M], but a more interesting approach is 
that of polariton physics [U El [6] , where the interaction with a nearby qubit can provide 
for both attractive and repulsive nonlinearities. 

The measurement of properties in this setup is also an interesting topic. By nature, 
the circuit-QED setups tend to be rather closed, as any additional cable or probe may 
be regarded as a source of decoherence. However, following Ref. [25], we envision an 
alternative that consists on a single transmission line running parallel to our setup of 
coupled cavities. Ensuring a weak coupling between both systems we can watch in real 
time the out-coupling of photons from the cavities into the line, thereby probing the 
frequency- and momentum-dependent correlations from Fig. [1} In addition to this, it 
would be interesting to probe the transport properties of these models, to see whether 
the pairing may assist in the mobility of photons, transmitting correlations through the 
chain, and probe the lack of insulating phases. 

5. Conclusions 

Summing up, in this work we have studied a lattice model of bosons coupled entirely 
through counter-rotating pairing terms. Unlike its relative the Bose-Hubbard model, the 
model from Eq. ([I]) does not experiment any quantum phase transition. Instead, for 
any value of the pairing, correlations are established among cavities, which we quantify 
using the pairing order parameter (aloj). We have shown that this pairing is maximal 
in momentum space, where it leads to the establishment of two-mode squeezing between 
different momenta. The selection of paired momenta can be controlled through the phase 
of the driving or the phase in Eq. ([I]). 

All these properties, including the numerically computed correlators and two-photon 
covariance matrix, can be explored in a setup that consists on superconducting resonators 
coupled by driven SQUIDs. An open problem which escapes the tools and methods shown 
in this work would be to study the transport properties of these models, to see whether 
the pairing propagates from the chain to the injected photons. 

Finally, we would like to point out that the model that we have studied could also 
be simulated using ultracold atoms in optical lattices. In this case the pairing term is 
provided by a Feshbach resonance that splits a molecule into its two fermionic or bosonic 
constituents. The split atoms would be trapped in an optical lattice which hinders their 
mobility and provides an on-site interaction U. The resulting model has exactly the form 



of Eq. (19), but demands a more complicated setup than the microwave photonics ideas 
shown above. 

JJG-R and AK acknowledge financial support from the European project 
PROMISCE, the Spanish MINECO project FIS2012-33022 and the CAM research 
consortium QUITEMAD S2009-ESP-1594. 
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Appendix A. Pertubation theory up to second order 

Motivated by the numerical result, we probe perturbatively the state 

N-l 

K) = (S)\^h (A.l) 

for an equally loaded lattice with uniform population n on each lattice site 
i = 0, . . . ,N — 1. Considering N lattice sites, the unperturbed hamiltonian, Hq, and 
the perturbative term, Hi, are given by 

Af-l N-l N-l 

Ho = uY^ 4ala,ai - /i ^ ajai , Hi = -J^ {e'^'44+i + H.c.) , (A.2) 

i=0 i=0 4=0 

with H = Ho + gHi. Its unperturbed energy is given by 

E^ = Nn [u{n-l)-fj] . (A.3) 

The state |u°) is not degenerate, hence we may refer to standard Schrodinger 
theory [26] for the perturbative ansatz. The corresponding expansion in state, 
\un) = K) + g \ui) + g^ \ul) + 0{g^), and energy, E^ = E^^ + g El + g^ El + 0{g^), 
provides no correction up to fist order in energy, E^ = 0. This is a direct consequence of 
the fact that Bloch states are now driven by excitations in pairs. The total number of 
excitations is no longer preserved. 

However, with 

we may obtain a correction for the eigenstate in first oder in g. Here, je^) refers to the 



rather trivial eigenstates of Hq in Eq. (A.2) with eigenenergies e^, i. e. Ho\es) = esjcs). 
Consequently the latter provides 



r N~l r N-l 

1=0 j=0 



with amplitudes 



77/ 

2U [n — 1) ~ jj, 
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and the notation |{ai, 02, ... , am.}i) marks the position i = 0, . . . ,N — 1 of the 
m-tupel ai,a2, ■ ■ ■ ,am in the otherwhise equally loaded state |m°). For ex- 
ample, taking into account periodic boundary conditions on the ring, it is 
\{n — l,n — IJat-i) = \n — l,n, . . . ,n,n — 1). 

From Eq. (A.5) a second order correction for the energy is immediately found to be 

El = j [nfn - (n + l)fn+i] . (A.7) 

Since there is no first order correction in energy, that is {u^\Hi\u^) = 0, one may 
rewrite the second order correction for the state as 

K) = y: ^i^i!f \^s) - 1 {f^ + fu K) . (A.8) 

Since we are interested in expectation values that lead to first and second order statistical 
moments, only the second term in the latter really matters up to second order in g. The 
first term requires at least higher moments in order to have an effect in the same order 
of g. We hence abdicate providing further details. 



